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ABSTRACT 



Aims. We study how IMF sampling affects the ionizing flux and emission line spectra of low mass stellar clusters. 
Methods. We performed 2 x 10 6 Monte Carlo simulations of zero-age solar-metallicity stellar clusters covering the 20 - 10 6 Mq mass 
range. We study the distribution of cluster stellar masses, M c i lls , ionizing fluxes, <2(H°), and effective temperatures, T^ s . We compute 
photoionization models that broadly describe the results of the simulations and compare them with photoionization grids. 
Results. Our main results are: (a) A large number of low mass clusters (80% for M c | us = 100 Mq) are unable to form an H n region, 
(b) There are a few overluminous stellar clusters that form H n regions. These overluminous clusters preserve statistically the mean 
value of (<2(H )) obtained by synthesis models, but the mean value cannot be used as a description of particular clusters, (c) The 
ionizing continuum of clusters with M dlls < 10 4 Mq is more accurately described by an individual star with self-consistent effective 
temperature (r* ff ) and <2(H°) than by the ensemble of stars (or a cluster T^ s ) produced by synthesis models, (d) Photoionization grids 
of stellar clusters can not be used to derive the global properties of low mass clusters. 

Conclusions. Although variations in the upper mass limit, m up , of the IMF would reproduce the effects of IMF sampling, we find that 
an ad hoc law that relates m up to M i us in the modelling of stellar clusters is useless, since: (a) it does not cover the whole range of 
possible cases, and (b) the modelling of stellar clusters with an IMF is motivated by the need to derive the global properties of the 
cluster: however, in clusters affected by sampling effects we have no access to global information of the cluster but only particular 
information about a few individual stars. 

Key words. Galaxies: star clusters, (ISM): Hn regions. 



1. Introduction 

H ii regions are gas clouds photoionized by the ionizing radi- 
ation of close-by stars and they are identified observationally 
with their characteristic emission-line spectrum. This spectrum 
depends on the physical conditions of the gas and the shape and 
intensity of the stellar radiation field, which depend in turn on 
the evolutionary status of the stellar component. However, in the 
general astrophysical context, we have access neither to the evo- 
lutionary status nor the physical conditions of the components, 
and we must derive them from observations. However, the over- 
all problem is difficult (and sometime impossible) to solve in an 
exact way because of the large number of unknowns involved, 
such as the exact 3-D gas and stellar distributions, the interre- 
lations and feedback between the different components, and in 
general, the unknown prior evolution of the system. 

One way to determine the physical conditions of the sys- 
tem is to perform a tailored analysis of the observed ob- 
ject (e.g.. ILuridiana et all I1999L iLuridiana & Peimbertl l200ll 
Luridiana et alj I2003I) . For large samples of data, it would be 
more economic (in terms of time and personal effort) to assume 
a set of hypotheses that simplify the problem as much as pos- 
sible but which, at the same time, allow us to obtain results as 
realistic and generic as possible. 

In terms of the physical conditions of the gas, a com- 
mon simplifying hypothesis is to assum e a given photoioniza- 
tion case, such as cases A, B, C, or D (Lurid iana et al.l 12009) 
for the observed object or set of objects. For the evolution- 



ary status of the stellar component, there are common sim- 
plifying hypotheses for the description of the ionizing con- 
tinuum produced by the cluster. Examples are methods using 
either individual st ars, such as the use of a cluster effective 
temperature (e.g., [S tasiriska 1990, Stasiriska & Scha e"rei1ll997l 
iBresolin et al. 1 119991 iBresolin & Kenicuttll2002h or a stellar en- 
semble where the different stellar components are mixed accord- 
ing to s ome specific prescriptions, as in evolutionary synthesis 
models (l01ofssonlll989l lG arcfa Vargas & Diaz 1994 are exam- 
ples of pioneering works and Dopita et al.ll2006l a more recent 
one). 

Although the ionizing continuum cluster may be more ac- 
curately represented by the ensemble of stars deterministically 
predicted by a sophisticated synthesis model, this may not neces- 
sarily be the correct approach. Evolutionary synthesis models as- 
sume average numbers of stars with different effective tempera- 
tures and luminosities determined by the stellar evolution and the 
stellar birth rate: the number of stars with a given mass that were 
born at a given time (parametrized by the initial mass function, 
IMF, and the star formation history, SFH). The averages cor- 
rectly describe the asymptotic propertie s of clusters with infinite 
(or a very large) number of stars ( see Cervino & Valls-Gabaudl 
120091 ICervino & Lu ridianall2005l [2006, for more detailed ex- 
planations). In more realistic cases, this proportionality is only 
valid as the average of several clusters, but it is not always a re- 
liable representation of individual clusters. In the case of low 
mass clusters (M c i us < 10 3 Mq), the number of stars is also 
not large enough to properly sample the IMF, which may pro- 
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Fig. 1. Connected dots: T* s vresus (vs.) stellar mass (m) relation 
for individual stars assumed in th is work. Open diamond s: T* s 
vs. stellar mass (m) relation from Diaz-Mill er et al.l (0*998) data. 
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Fig. 2. r* ff vs. Q(He°) / Q(H°) ratio for individual stars assumed 
in this work. 



duce biased results when the observations are analysed by com- 
paring them to the mean value obtained by synthesis models 
dCervino & Valls-Ga baud 2003). These effects are expected to 
exist in a great part of the ionizing clusters, since the major- 
ity of the embedded clusters have low masses (M c i us < 10 3 Mq, 
lLada & Ladall2003l) . 

Given the decay of the IMF at high masses, the net effect 
is that a large number of low mass clusters (but not all) will 
not produce massive stars. In these cases, we demostrate that 
the spectrum of an individual star can provide a closer fit to the 



cluster ionizing continuum than the result of a synthesis model 
used in a deterministic way. 

In this paper, we show, by means of Monte Carlo simulations 
of zero age stellar clusters, that between these two approaches 
(individual star versus (vs.) ensemble representations) there is a 
smooth transition that depends on the cluster mass. We establish 
how representative individual stellar spectra are of low cluster 
masses and show their implications for estimating the physical 
properties of clusters by means of their emission line spectra. 

In section 2, we present our Monte Carlo simulations and 
probability computations. In section 3, we describe the photoion- 
ization simulations and analyse the results by means of diagnos- 
tic diagrams. The implications of the analysis of the results is 
discussed in section 4. Finally, in section 5 we summarize the 
main conclusions of this work. 

2. The distribution of the ionizing properties of 
clusters 

The first step in characterizing the emission line spectra of H n 
regions is to attempt to reproduce their ionizing continuum. For 
this purpose, we used the cluster effective temperature, T c ™, 
instead of a detailed spectral energy distribution. The cluster 
is defined as the effective temperature of a reference star 

(T*J ef ) that has the same ratio g(He°)/G(H°) as the cluster 
(Mas-Hesse & Kunthlll991h . i.e., the same number of ionizing 
photons emitted below the He i Lyman limit at 504 A, Q(He°), 
with respect to the total ionizing photon number, Q(H°). 
The reference scale was defined in the following way: 

1. The relation between Q(H°) and initial mass was obtained 
using the ZAMS values published by IDiaz-Miller et al.l 
( 1998) (their Table 1), after excluding their two lowest mass 
points (see below). 

2. Low mass stars (m < 2OM0): 

(a) 

provided bv lDfaz-Miller et all dl998h 
(2(He°) values were derived from t he corresponding 
Q(H°) , \ogg, and r* ff values from IDiaz-Miller et al.l 
( 1998) assumin g the conti nuum shape given by ATLAS 
atmospheres (Kurucz 1993). 
(c) A 0. 15 M0 star was included to cover the complete range 
of stellar masses implicit in the IMF. Th e T*g B is the 
r* ff from a 0.15 M0 star in the tracks by iGirardi et al.l 
(2000). The g(H°) and g(He°) values for this star were 
obtained by linear extrapolation in the log m - log Q(X) 
plane of the points for 2.6 and 2.9 Mq obtained in (a) and 
(b) above . The reason for not using the two lowest mass 
points of IDiaz-Miller et al.l dl998h is the poorer agree- 
ment of the linear extrapolation in the log m - log T* ff 
plane from those points with the r* ff of a 0.15 M0 star. 

This can be seen in Fig. Q] w here the m — T** e f rela- 
tion used in this work and in Diaz-Miller et al. d!998h 
are shown. 
High mass stars (m > 20Mq): 

(a) r *' ref values w ere obtained directly from Cloudy's 
(Ferland et al. 1998) implementation of CoStar 



values were assumed to be the T* ff for ZAMS stars 



(b) 



mp li 

Schaerer & de Koterj 1 19971) atmospheres, and the 



associated evolutionary tracks. The Cloudy inputs were 
the initial mass and the age (in our case 0.05 Ma for a 
ZAMS representation) and provided the corresponding 
atmosphere model and T* s from a combination of the 
CoStar grid. 
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Fig. 3. The color-coded area shows the M c i us 

the plot is the m - T*^ ef reference scale plotted in Fig.Q] Note that the stellar m - T*^ 1 reference scale defines the low mass limit 
of the M c i us - T^ s distribution. The solid-line contour on the right of the plot, reaching M c i us higher than 10 4 Mq, corresponds 
to the simulations of the high-mass cluster set (c.f. Fig. @). The inner box shows a detail of the region around 100 Mq, where the 
maximum value of the cluster T^ s is reached. 



(b) The Q(He ) values were derived from the resulting at- 
mosphere model assuming a continuum level given by 
the corresponding log g(H°) obtained in step 1 above. 



The resulting m - T*J ef relation and r* ff ref vs. £(He°) / Q(H°) 
reference scale are shown in Figs.Q]and|2j respectively. 

We note that the r*^ vs. g(He°) / Q(U°) r eference scale 
depends on the stellar atmosphere code (see iMartins et al.l 
l200l ISimon-Dfaz & Stasihsk^Hool for an illustrative study). 
Although this choice is critical to the study of the detailed 
emission-line spectra of individual stars and stellar clusters, it 
has a minor impact on the present study. Our goal is to estab- 
lish a reference scale to allow a self-consistent comparison of 
the ionizing properties of individual stars and clusters. Any ref- 
erence scale is valid as long as it provides a non-degenerate re- 
lation between the individual stellar mass (and age) and Q(He°), 
(2(H°), and their ratio (and, hence, T*^ et ). With this is mind, we 
used ATLAS and CoStar model atmosphere grids,partly also be- 
cause they are directly available in Cloudy and are the ones used 
by the evolutionary synthesis code considered in our study. 



2.1. Method 

We computed two sets of simulations. The first set was com- 
posed of 10 6 Monte Carlo simulations of clusters with masses 
from 20 to 10 4 Mq (in the following, the low-mass cluster set). 
The second set consists of 10 6 Monte Carlo simulations of clus- 
ters with masses from 10 3 to 10 6 Mq (high mass cluster set). We 
note that there is an overlap between the mass range of the two 
sets of simulations. 

Preliminary cluster masses, M c \-o, were obtained by means 
of random sampling of the initial cluster mass function (ICMF) 
in the cluster mass ranges defined for each set. The assumed 
ICMF has the form dN/dM c \- oc M d 2 for all the mass clus- 
ter range considered. This rep resentation of the ICMF was pro- 
posed bv lLada & Ladal (2003) for clusters with masses between 
~50 and 10 3 Mq. The application of this ICMF t o the complete 
mass r ange used here is suppo rted by the studies of Zhang & Falll 
(I1999I) and lHunter et al.l (I2003I) . 

Once a cluster mass was assigned, we used a Salpeterl(fl95l 
stellar initial mass function (IMF) with a lower mass limit mi ow = 
0.15 Mq and an upper limit of m up = 100 Mq, to determine the 
masses of the individual stars, m, in the cluster by performing 
random sampling of the IMF. The sampling of the IMF continues 
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Fig. 4. The color-coded area shows the M c i us - T^ s distribution for the high-mass cluster set. The solid-line contour on the left, 
reaching M c i us lower than 10 3 Mq, corresponds to the low-mass cluster set (c.f. Fig. [3]). 



until the preliminary cluster mass M c i ; o is reached. Hence, all 
simulated cluster have masses M c i us higher than the preliminary 
Mcijo and the corresponding ICMF of a given set of Monte Carlo 
simulations does not follow exactly a power law with exponent 
-2. 

In each individual cluster, the Q(H°) and Q(He°) values of 
each star of mass m were obtained by linear interpolation of the 
log m - log Q(X) relation obtained previously. We obtained the 
cluster's integrated Q(H°) and Q(He°) by adding the contribu- 
tion of each individual star, and obtained the cluster's T™^ from 
a linear interpolation in the log g(He Q )/<2(H°) - relation 
shown in Fig. [2] 

We rebin neither the stellar masses nor the cluster masses, so 
our results cor rectly map the underlying d istribution of physical 
properties (see lCervino & Luridianal2004 for different ways of 
implementing the Monte Carlo simulation). 

2.2. Results 

We now present the resulting multi-dimensional distribution fo- 
cusing on M c i us , the number of stars in the cluster (N c \ us ), T°^, 
and Q(H°). We compare the results of both sets of simulations 
using, when possible, the cluster distributions with the reference 
scale assumed. We also analyse the statistical properties of the 
distribution of Q(H°) and r c I us as a function of the cluster mass. 



2.2.1 . The M clus - Tf» s distribution 

The color-coded area in Fig. [3] shows the bi-parametric M c ] us - 
r^ us distribution of the low-mass cluster set. For comparison, the 
contour of the M c i us - Tf^ s distribution of the high-mass cluster 
set is overplotted as a black solid line. Figure |4] shows the same 
distribution for the high-mass cluster set but with the opposite 
coding (color-coded for the high-mass cluster set, black solid 
line for the low-mass cluster set). 

The sharp vertical edge in the distributions is due to the min- 
imum cluster mass assumed in the simulations sets (20 M0 for 
the low-mass, and 10 3 Mq for the high-mass cluster sets). The 
solid line with dots on the left-hand side of Fig.[3]is the m - T** et 
reference scale, already plotted in Fig.[T] 

The m — T*£ ef reference scale coincides with the low mass 

ett 

envelope of the M c i us - r^ us distribution for M c i us lower than the 
upper limit of the stellar IMF (ra up ). This result is obvious since 
the M c ius - 7^f US distribution is limited by the extreme cases of 
clusters composed by a individual star with mass m = M c i us (see 
also Fig. [5] below). The corresponding T^ s of a cluster com- 
posed of a mixture of stars is not so obvious. In general, the re- 
sulting r^ us of a cluster is between the minimum and maximum 

r*^ ef of the stars present in the cluster. However, its exact value 
depends on the relative luminosity of the cluster components. A 
cluster containing only a 50 Mq star will have a T^ s similar to a 
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Fig. 5. The color-coded area represents the (N c \ m , Tf^ s ) distribution of the low-mass cluster set. The solid-line contour, extending 
to A^ c i us larger than 10 5 , shows the distribution of the high-mass cluster set. The vertical discontinuity that appears at N c i us = 10 2 is 
a visual artifact caused by a change from linear binning to logarithmic binning. 



cluster with a 50 M0 star in addition to about 20 stars with 5 Mq 
(M c i us = 150 Mq), since the g(He°) and Q(H°) produced by such 
a cluster is dominated by the 50 M0 star. However, a cluster with 
a 50 Mq star will have a different to a cluster with both a 
50 Mq star and a 100 Mq star since both have a non negligible 
contribution to Q(He°) and Q(H°) (although this particular com- 
bination has a low probability of existing given the IMF). The 
situation in the case of the maximum r^ us (= r*^ ef ) is illustrated 
in the inner box of Fig. [3] The figure shows the locus of individ- 
ual simulations in the region around the maximum value of the 
cluster T^ s . The dashed line shows the value of the maximum 

r*^ ref , equal to 5.21 x 10 4 K. As expected, there are no simula- 
tions in which I™ 8 is higher than this value but, for the assumed 
ICMF, there is a sizeable number of clusters in the 100 - 10 3 Mq 
range with such a high effective temperature, i.e. with at least a 
100 Mq star that dominates the ionizing flux. As the relative 
contributions of the stellar mixture becomes more homogeneous 
among clusters, the T^ s distribution becomes increasingly nar- 
row, and eventually collapses to t he asymptotical value T^ s 
obtained by synthesis models (see ICervino & Luridianal 2006, 
ICervifio & Valls-Gabaudl2"009llBarker" et al.ll2008l for more de- 
tails), which corresponds to 4.54 x 10 4 K for the assumed age, 
IMF, and m - T*i ef reference scale. 



The m - T*^ e reference scale not only defines the low mass 
envelope of the M c i us -r^ s distribution, but also produces multi- 
modality in the distribution. Steeper slopes in the m - T*^ et refer- 
ence scale correspond to less populated zones in the M c i us - T^ s 
distribution. As an example, the steeper slope in the 17.5 - 20 
Mq range of the m - T*£ ef reference scale, corresponding to the 
3 - 3.5 x 10 4 K range, is reflected in a minimum in the density 
of cluster simulations. 

2.2.2. The N clus - T^ s distribution 

Figure|5]shows the N c i us - T^ s distribution for both sets of sim- 
ulations (only the contour of the simulations covered by the the 
high mass cluster set is shown). Both set of simulations produce 
a pyramidal distribution shape. 

The left part of this pyramid is produced by the limitation 
imposed by the minimum cluster mass in both sets. The clusters 
that define the lower N c i us envelope in the low and high mass 
simulation sets have masses of between 20 and 100 M (low 
mass set) and 1000 M (high mass set). The low A^ c i us end of the 
simulations shows that the clusters with only a few stars have in 
fact a high Tf^ s : they are formed by only a few high mass stars 
due to the boundary conditions imposed on the simulations. 
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Fig. 6. Color-coded distribution of the cluster ionizing flux (Q(H )) and T^ s for the low- and high-mass cluster set. The shaded 
contour represents the low-mass cluster set, the solid-line contour the high-mass cluster set. A strong correlation between both 
quantities is seen for the low-mass cluster set case. The same corelation can be seen in the lower part of the distribution of the 
high mass set, and is lost in the upper part, where the distribution tends to the asymptotic T^} s . The insets show the histograms of 
the 7^ ls distribution for two values of (2(H°). The range in Q{H°) for those histograms are marked with vertical dotted lines. As a 



reference, we also plot the 2(H ) 



r*^ ef reference scale (line with symbols). 



The right part of the pyramid is controlled by the asymp- 
totic trend of the clusters with larger N c \ us towards the asymp- 
totic r c Jj ls value of clusters with infinite number of stars. 



2.2.3. The g(H°) - Tf™ distribution 

In Fig. [6] we show the Q(H°) - 7^ us distribution for both sets 
of simulations, with the same coding as adopted in Fig. [5] As a 
reference, we have also plotted the Q(H°) - T*'^ e relation (black 
solid line with dots), which defines the upper envelope of the 
distribution in the low-mass set of simulations. The 2(H°)-7'^ US 
relation is very narrow for simulations with M c i us < 10 4 Mq. 
Moreover, for the low-mass set given and a fixed Q(H°), the most 
probable T~j£* value is the one corresponding to an individual 
(reference) star with this given Q(H°), i.e., the most massive star 
compatible with this Q(H°) value. This is shown by the small 
histograms of T°^ s for simulations with log Q(H°) around 44 (as 
an example of low Q(H°)). 

At intermediate log Q(H°) values (for example 48.58, the 
value from synthesis models for a 100 Mq cluster also shown in 



the figure) there appears to be both bimodal distributions and a 
larger scatteiQin the Q(B°)-T^ S distribution. At such log Q(H°) 
values, there are clusters containing either an individual domi- 
nant star or a mix of stars that mimic clusters with a well sam- 
pled IMF but with a lower m up value (i.e. high-mass-star defi- 
cient clusters). 

Finally, when the Q(H°) of the cluster is larger than the 
(2(H°) of any individual possible star (in this case, clusters with 
masses higher than 10 4 Mq) the cluster T^ ls distribution slowly 
converges to the asymptotic T^ s value, obtained by synthesis 
models that do not take into account sampling effects. 

The figure illustrates the effects of sampling on observable 
quantities. In the low Q(H°) regime, if we know the Q(H°) of a 
low mass cluster, we can assign a 7^ us with some confidence, 
since both relate to similar (dominant) stars, but we are unable 
to determine the global properties of the cluster (M c i us or A^ c i us ) 
with similar confidence. We can, however, determine informa- 
tion about particular stars in the cluster (the most luminous 

1 The large area covered by the distribution in this range indicates that 
there are two relatively narrow distributions centred on different Tfg s 
values, rather than a single unimodal distribution with a larger scatter. 
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ones). Th is is the case for cluster s below the lowest luminos- 
ity limit (ICervino & Luridiana 20040, where extreme sampling 
effects occur, information about only these particular stars is ac- 
cessible, not for the cluster as a whole. 

Near the lowest luminosity limit, there are clusters dominated 
by individual stars (defining a narrow distribution peaking near 
the assumed reference scale) and there are clusters whose emis- 
sion is dominated by an ensemble of stars. This second kind of 
clusters can be characterised as clusters deficient in high mass 
stars. In this case, sampling effects can fool us when the asymp- 
totic values obtained by synthesis models are used to obtain 
global properties (M c i us , N c \ w , or ages). When we modify m up 
in the IMF to mimic the low luminosity clusters, we lose a size- 
able number of clusters dominating an individual massive star. If 
we consider clusters dominated by an individual massive star, we 
lose a sizeable number of clusters deficient in high mass stars. 

In the regime of large Q(H°) (values larger than around 10 
times the lowest luminosity limit) both T^ s and Q(H°) are de- 
fined by the proportionality of stellar types provided by the IMF 
with the correct m up , so they provide information about the clus- 
ter as a whole and variations in the numbers of particular stars 
have little impact on the integrated luminosity. 

This strong correlation is present only for observed quan- 
tities and is lost when the cluster mass is used: the Q(H°) vs. 
M c i us distribution (not plotted) shows a similar scatter as the 
M c i us - T^? s distribution shown in Figs.[3]and[4] 



2.2.4. Statistical properties of Q(H°) 

Figure [7] shows the distribution of log Q(H°) in clusters with 
masses around 100 Mq. The position of the mean value of the 
distribution ((Q(B°)} = 1.1 x 10 48 ph s _1 , log<0(H°)> = 48.04) 
is plotted as a vertical dashed line. We also determined the 
mean Q(H°) for other mass ranges, which consistently pro- 
duces the same number of ionizing photons normalized to mass 
<2(H°)> = 1.9 x 10 46 ph s " ] 

yields <fi(H )) = 1.9 X 10 48 ph s 
discrepancy between this value and the one obtained from the 
distribution around M c i us = 100 Mq is consistent with low mass 
clusters being more numerous because of the ICMF 

Approximately 84% of the clusters have values below the 
mean, and only 16% above it. The properties of these overlu- 
minous clusters, and therefore their contribution, depend on the 
upper limit mass of the IMF, m up , because a change in the upper 
mass leads to a change in the underlying probability distribu- 
tion (the very IMF). Therefore, the average value of quantities 
such as (Q(H°)) will be not preserved during changes in m up . In 
this sense, although there is a strong (statistical) correlation be- 
tween M c i us and the star of maximum mass m max in the cluster, 
and an average (mm M ) c an be obtained as a function of M c i us 
dWeidner & Kr oupa 2006), it is incorrect to use (m m3X (M c \ us )) as 
an upper mass limit (m up ) when defi ning the IMF (whose up- 
per mass limit m» v appears constant: IWeidner & Kroupal l2004. 
Mafz Apellaniz et al.l l2007T) . Not only does this change fail to 
preserve (Q(H ))): it fails to preserve even (m max (M c i us )). 

This experiment allows us to illustrate some issues in the 
modelling of stellar clusters: once an observable is fixed (the to- 
tal mass of the cluster in this case), the other observables (the 
number of stars, Q(H°), or T~£* in our case) vary following 
given probability distributions. In other words, typical relations 



M- 1 (log<2(H u )> = 46.27), and 
1 for M clus = 100 M Q . The 



2 The lowest luminosity limit is defined as the luminosity of the 
brightest individual star that can be present in the cluster for the given 
evolutionary conditions. 
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Fig. 7. Distribution of log g(H°) for simulations with 90 Mq < 
M c ius ^ HOMq obtained from the low mass cluster set. The 
probabilities were computed using the simulations in the given 
mass range (40 502 simulations) and found to have a small 
dependence on the ICMF. The vertical line shows the mean 
<G(H°)> value (log<g(H°)> = 48.04 of the distribution. 



obtained by synthesis models should be considered valid as a 
mean of the distribution of observed clusters, but not necessarily 
valid for particular ones. In addition, there are strong correla- 
tions between different observables (such as Q(H°) and 7^ ls ), 
which reduce the scatter in the distributions considerably, espe- 
cially for undersampled clusters. As a drawback of these situ- 
ations, the observables do not reflect the cluster properties, but 
just the properties of the most luminous stars in the cluster. 

2.2.5. Probabilities for given mass ranges 

We now consider the distribution of the T^ s for some spe- 
cific mass ranges. Figures [8] [9] and [10] show three cuts of the 
distribution for total cluster mass values around 100 Mq, 10 3 
Mq, and 5 x 10 5 Mq, respectively. The masses considered are 
90 Mq < Mcius < 1 10 Mq for Fig. [8] which cover 40 502 sim- 
ulations. In the case of Fig. [9] which includes simulations from 
the two sets, we have used 900 Mq < M c i us < 1.1 x 10 3 Mq in 
the low mass cluster set with 4 091 simulations, and 10 3 Mq < 
M clus < 1.1 x 10 3 Mq in the high mass cluster set with 90621 
simulations. Finally, the range in Fig. [TO] is 4.5 x 10 5 Mq < 
M c j uii < 5.5 x 10 5 Mq with 373 simulations. In all cases, a renor- 
malization to the number of simulations in the considered mass 
range was performed for an easy comparison. 

The T^ ls distribution shown in Fig. [8] span a wide range of 

values corresponding roughly to the T*'™ 1 of stars with masses 
between 2 and 100 Mq. The distribution is bimodal with a mean 
value of around 2 x 10 4 K, which is clearly biased with 
respect to the asymptotic value of 4.54 x 10 4 K obtained from 
the average of the Tfg S of clu sters with M c i us > 10 5 Mq (see 
ICervino & Valls-Ga baud 2003, for more details of this kind of 
bias where the observable is obtained from a ratio). 
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Table 1. Probabilities for the simulations with M as 100 Mq. 



interval (10 4 K) 


# Simulations 


Probability (%) 


Partial probability (%) 


Stellar mass (Mq) 


4.92-5.51 


259 


0.64 


3.31 


100 


4.24 - 4.92 


2410 


5.95 


30.80 


50 


3.70-4.24 


3503 


8.65 


44.77 


25 


3.4-3.7 


1653 


4.08 


21.12 


20 


0.2-3.4 


32677 


80.68 




< 20 



Notes. The first column shows the r^ us interval in which probabilities are computed; these ranges are also shown in Fig.[8] The second column gives 
the number of simulations in each T?^ 8 range. The third column shows the probability of each T^ s range when the complete set of simulations 
is considered. The fourth column shows the corresponding probability when only the subset of simulations that can generate an Hn region is 
considered. The last column lists the stellar mass corresponding to the middle point of each temperature range. 




T^(K) 

Fig. 8. Distribution of T^ s for simulations with 90 M Q < 
M c i us < 1 10 Mq obtained for the low mass cluster set. The prob- 
abilities are relative to the number of simulations in the given 
mass range (40 502 simulations). Vertical lines delimit the re- 
gions used for the computation of the probabilities in Table [TJ 



For M c ius around 10 3 Mq, the shape of the distributions is 
also bimodal, but the mean is more consistent with the asymp- 
totic value. As expected, this bimodality disappears when M c i us 
is a round 5 x 10 5 M(?,, where only a small asymmetry remain s 
(cf . ICervino & Valls-Gabaudll2003i ICervino & Luridianall2006l) . 

These histograms allow us to determine the clearest rep- 
resentation of the ionizing continuum by performing a non- 
exhaustive exploration of the T^ s distributions. In the follow- 
ing, we describe how focus on the histogram with M c i us around 
100 Mq, which corresponds to the most n umerous stellar clus - 
ters according to the ICMF proposed bv lLada & Ladal (|2003). 
We divided the 7™ s range into five intervals. The middle point 
of the first four intervals correspond roughly to the T* ff of ZAMS 
stars with 100, 50, 25, and 20 Mq, the last one including all 
the clusters with T c ^ s lower than 3.4 x 10 4 K, which we use as 

ett ' 

an approximate lower limit for the developing of an observable 
H ii region. In Table Q] the number of simulations and the corre- 
sponding probabilities for each interval are shown. 

Table[T]shows that most of the zero-aged clusters of 100 Mq 
(around 8 1 %) are unable to generate an H n region because they 




3-10 4 4«10 4 5»10 4 

T?»* (K) 

Fig. 9. Distribution of T^ us for simulations with 900 M < 
M c j us ^ 1.1 x 10 3 Mq obtained from the low-mass cluster set 
(4 091 simulations) and 10 3 M o < M clus < 1.1 x 10 3 M o , ob- 
tained from the high-mass cluster set (90 621 simulations). 

do not have sufficiently massive stars. In the remaining clusters, 
around 66% (13% of the total) have T*Q* below the asymptotic 
value of 4.54 x 10 4 K and only 34% (6% of the total) have a 
value comparable to or larger than the asymptotic one. 

3. Effects in H n regions 

The dispersion in T^ s shown in the previous section implies that 
there is also a dispersion in the cluster energy distributions. We 
now show in detail how this dispersion affects the emission-line 
spectra produced by low mass stellar clusters. We also show how 
it can render useless the predictions of photoionization grids of 
clusters modeled as ensembles of stars such as those produced 
using synthesis models in a deterministic way, e.g., a univocal 
value of the ionizing flux for a given cluster mass. 

We focus on low mass clusters with M c i us around 100 Mq, 
which are the most abundant ones. We performed photoioniza- 
tion simulations of this kind of clusters in a simplified way by 
considering: clusters with 5 ionizing stars (N*j on ) and the ioniz- 
ing continuum of a 20 Mq ZAMS stars, clusters with N*j on = 4 
and the ionizing continuum of a 25 Mq ZAMS star, clusters with 
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Table 2. Photoionization simulations for clusters with M c i us 
100M©. 







Stellar mass 




Complete IMF 




100 


50 


25 


20 


(0.15-100) 


N, 


l 


2 


4 


5 


192.6 


N *,ion 


l 


2 


4 


5 


0.23 


log(fi(H°) 


49.87 


49.34 


48.84 


48.58 


48.58 


r^io 4 K) 


5.01 


4.53 


3.81 


3.53 


4.60 


Table 3. Photoionization simulations with log 


Q(H°) = 49.87. 






Stellar mass 




Complete IMF 




100 


50 


25 


20 


(0.15-100) 


N, 


1 


6.76 


42.66 


97.72 


1778.72 


N ■ 


1 


6.76 


42.66 


97.72 


2.1 


log(G(H°) 


49.87 


49.87 


49.87 


49.87 


49.87 


r£' s (io 4 K) 


5.01 


4.53 


3.81 


3.53 


4.60 



Notes, log Q(H°) equals to 49.87 corresponds to the Q(H°) of an indi- 
vidual 100 M© star. 



Teff S (K) 



Fig. 10. Distribution of 7™? for simulations with 4.5 x 10 5 M Q < 
A^cius ^ 5.5 x 10 5 Mq obtained from the high-mass cluster set 
(373 simulations). 




Energy (eV) 

Fig. 11. Ionizing continua for the scenario of fixed log Q(H°) = 
49.87. Different lines correspond to different shapes of the ion- 
izing continuum, as explained in the text. Note the small differ- 
ences among the results of synthesis models (labelled as IMF) 
and the continuum of a 50 M0 and 100 Mq star, consistent with 
the similarity in between the cluster and star. 



N*,ion = 2 and the ionizing continuum of a 50 M0 ZAMS star, 
and clusters with N % i oa = 1 and the ionizing continuum of a 
100 M0 ZAMS star. All ionizing continua have been obtained 
from CoStar models (ISchaerer & de Koterl ll997). which consis- 
tent with the in - r* ff relation described in Sect. [2] The summary 
of this set of photoionization simulations is presented in Table 



[2] We recall that this situation covers only 20% of all possible 
cases, since 80% of clusters of this mass do not produce stars 
massive enough to develop an H n region (cf. Table[TJi. 

As a reference value, we also used the ionizing contin- 
uum of a solar-metallicity zero-a ge stellar population synthe- 
sis model compute d with SED@ dCervino & Mas-Hessd 1994, 
ICervino et al . 2002E which uses stellar atmospheres identical to 
those adopted by ourselves in Sect. [2] However, SED@ assumes 
a m - 2(H°) relation given by evolutionary tracks, which is dif- 
ferent from that used in this work, so that there are some small 
discrepancies. The most important difference is that SED@ as- 
sumes m up = 120 Mq, hence slightly larger asymptotic val- 
ues are expected in both Q(H°) and T^ s . SED@ provides a 
T^ s = 4.59 x 10 4 K, which is 500 K higher than the asymp- 
totic value obtained by the simulations. SED@ also infers that 
log<e(H°)> = 3.8 x 10 46 erg s" 1 M" 1 , which also differs from 

the asymptotic value of 1.9 x 10 46 erg s _1 Mq 1 obtained here. 
However, these discrepancies can be explained by the different 
value adopted for m up in the simulations and SED@, and do not 
affect the current discussion. 

When synthesis models are used, it is necessary to distin- 
guish between the number of stars in the cluster, N*, and the 
number of ionizing stars in the cluster, N t j on . These numbers 
can be obtained by simple integrals and do not depend on the 
synthesis model. The mean number of N,j on for a 100 M0 clus- 
ter is 0.23, consistent with the result of Table [TJ where ~ 20% of 
the clusters have stars able to develop an H 11 region. 

We performed additional photoionization simulations for a 
fixed Q(H°) = 49.87, corresponding to the Q(H°) of an individ- 
ual 100 Mq star. The summary of these sets is shown in Table 
® 

As we have already discussed, the tables again show that, in 
a realistic case, once the cluster mass is fixed, the other prop- 
erties of the cluster (N t , N*; lon , and Q{H°) as examples) have 
distributed values. In the case of M c i us = 100 Mq, there is no so- 
lution that would simultaneously satisfy the results in continuum 
shape, continuum intensity (parametrized as Q(H )), and cluster 



3 The SED@ population synthesis models are available at 
http://www.laeff.inta.es/users/mcs/SED with values renor- 
malized to the used 0.15-100 Mq IMF mass range. 
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mass/number of stars described by the average value obtained 
from synthesis models. Given that the continuum shape and lu- 
mino sity control the emiss ion line intensity in a non-linear way 
(see IVillaverde et alJl2009l for more details), using an average 
ionizing continuum of clusters in function of the cluster mass 
as input of photoionization codes, as synthesis models provide, 
gives unrealistic results. 

The models with fixed Q(H°) allow us to compare the ioniz- 
ing flux of the different cases within a reference scale. FigureQj] 
shows the ionizing continuum for all the ionizing fluxes calcu- 
lated for the Q(H°) fixed case. The shape of the average continua 
obtained with synthesis models is similar to that of clusters with 
only 50 Mq or 100 Mq stars : the asymptotic cluster 7^ us is 
between the r* ff of a 50 M Q and that of a 100 M Q ZAMS star. 
The figure also shows the large differences above 24.6 eV in the 
average continuum obtained from synthesis models and the sim- 
ulations with 20 and 25 Mq stars. 

The set shown in Table [3] shows us another interesting con- 
clusion. Synthesis models require about 1800 stars (i.e. cluster 
masses around 10 3 Mq) to produce Q(H°) similar to an individ- 
ual 100 M0 star. This allows us to establish a mass limit for the 
use of synthesis models in stellar clusters: clusters with masses 
lower than 10 3 M0 are less luminous than the individual stars 
that the cluster model contains when described as the average 
obtained by synthesis models (the lowest luminosity limit es- 
tablished bv lCervino & Luridi ana 2004, lCervino et al.l2003l) . As 
we have shown (c.f. Fig. [6] and discussion in the text), the low- 
est luminosity limit indeed establishes only the minimum cluster 
mass value to avoid bolometric inconsistencies, although impor- 
tant sampling effects are found for mas ses 10 times higher th an 
the lowest luminosity limit value dCervino & Lurid iana 2004). 

3.1. Effects on the emission line spectrum 

We now study the effects of using an ionizing continuum pro- 
duced by a cluster consisting of individual stars instead of the 
average spectrum obtained by synthesis models. We inputted 
the different contin ua into the photoionization code Cloudy 
dFerland et al.l 1 19981 version C08). In addition, we assumed 
a rather simple description of the gaseous component: an ex- 
panding radiation-limited nebula with a 1 pc inner radio, density 
equal to 100 c m" 3 , and the same solar metallicity mixture as in 
Table 1 from Dopita et al. (2006), although without taking into 
account any depletion of metals. 

To analyse easily the results obtained , we used diagnosti c 
diagrams and compared with the grids by Dopita et al.1 (12006b . 
These grids were parametrized as a function of a parameter %, 
which is the logarithm of the ratio of the cluster mass to the am- 
bient pressure. For comparison, we chose the % that most closely 
reproduces the position of our Cloudy simulations using the av- 
erage ionizing continuum obtained by solar metallicity zero-age 
synthesis models (% = -2 for models with fixed M c i us , and H = 2 
for models with fixed g(H°) ). 

Figure [T2l shows the class ical diagram log([Q inj/l 5007 /H/j) 
versus log([N n]A 6584/Ha) ([Baldwin et al.ll 19811) for the two % 
values considered. The figure shows large discrepancies when 
the predictions of the grids are compared with the results of clus- 
ters formed with particular stars. In particular, the models with 
ionizing continua described by 100 and 50 Mq stars and fixed 
M c i us are outside the grid coverage in a region that would cor- 
respond to ages younger than those described in the grid (cf. 
Fig[T2]left). The situation is better for clusters with fixed Q(H) 
(cf. Fig[T2]right) where the results of synthesis models (labeled 
as IMF) coincide with the results of clusters with only 50 M0 



stars. In both cases, clusters with ionizing continuum described 
by 25 and 20 M0 stars are in the region covered by the grid, 
although, if the grid is used to estimate ages/metallicities, they 
are confused with old clusters with a lower metallicity than the 
one they actually have. The error in the age estimation can be 
easily explained: the average ionizing continuum produced by 
synthesis models has a large contribution from the stars at the 
main-sequence turn-off (the hotter ones in the cluster). As time 
evolves, the stars at the turn-off become less massive and cooler, 
so the average cluster ionizing continuum is more accurately de- 
scribed by the 25 and 20 Mq stars. For clusters affected by sam- 
pling effects, the cluster is not actually older, but just deficient 
in massive stars. This effect would also be reproduced by arti- 
ficially imposing a lower limit of the m up of about 20 M0: in 
these cases, the integrated light of the cluster will not vary be- 
tween zero age and the age when stars with m up begin to leave 
the main-sequence, and a zero age is compatible with the model 
result. However, although it would seem to be a good solution 
for some low-mass clusters (e.g., of 100 Mq), this variation in 
m up is not valid in general (7% of the clusters have at least a star 
more massive than 50 Mq, and 80% of clusters have a maximum 
m in the particular realization of the IMF that has no stars able 
to produce an H 11 region). 

n — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — J — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — J — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — J — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ri 



R=-2 




-2-10 1 
Log([NII]A,65847 [OH]A3727) 



Fig. 13. Log([Om] 5007 A//[On] 3727, 3729 A) versus 
Log([Nn] 6584 A/[On] 3727, 3729 A) for all the simulations. 
The figure only shows the grid for 7? = -2, which corresponds 
to the cases with fixed M c i us . Lines and symbols as in Fig. Q~2] 
left. Note that the model of 100 Mq for fixed M c i us and the cases 
of synthesis models (labeled as IMF), 50 and 100 M0 for fixed 
Q(H°) are in the same position on the plot. 

A similar situation can be seen in Fig. [T3l 
which shows log([Ora]^4959,5007/[Oii]^3727) vs. 
log([NiiU6 584/rOnl/l3727). This diagram was firstly 
proposed by Dopi ta et al.l d2000) for metal licity and ionization 
parameter estimation and was also used by Dop ita et~ai] (2006) 
for age estimations. The figure only shows the grid for fi = -2, 
which corresponds to the cases with fixed M c i us . The vertical 
axis again shows a large dispersion because of the differences 
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Fig. 12. log([Om] 5007 A / H/3) vs. log(|"N nj 6584 A/ Ha) dia gram for the individual-star zero-age solar metallicity simulations. 
The results are compared with the grid by iDopita et al.1 (2006) where dotted lines join points of equal ages, and filled lines join 
points with equal metallicity. Left: Grid with logarithm of the ratio of the cluster mass to the ambient presure, %, equals to -2, 
simulations with fixed M c i us = 100 Mq are plotted as large filled symbols (and models with fixed Q(H°) with small open ones). 
Right: Grid with % — 2, simulations with fixed log Q(H°) = 49.87 are plotted as large open symbols (and models with fixed M c i us 
with filled open ones). 



in the ionizing continuum of the simulations, which can be 
confused with age variations. On the other hand, the dispersion 
on the horizontal axis is very small, indicating that the diagram 
really depends on the gas metallicity. 

Variations in the 'R parameter are expected that for each of 
the particular cases to porvide a better result. By the same argu- 
ment, an even better result would be obtained by performing a 
customized analysis for each cluster!, taht would only provide 
information about particular stars in the cluster, but not about 
the global cluster properties, as the grids do (excluding the cases 
of rough metallicity estimations). In conclusion, grids of pho- 
toionization models are not useful enough tools for deriving the 
properties of stellar clusters affected by sampling effects. 

4. Discussion 

We have desmostrated the limitation of the use of average val- 
ues obtained by synthesis models to describe the properties of 
individual clusters with a direct application to H n regions. Of 
course, this limitation of the average is proportionally reduced as 
the number of stars in the observed cluster increases. However, 
the situation fo r low mass stellar clusters (the most abundant 
ones, following lLada & L ada 2003) is complex and the use of 
the mean value obtained by synthesis models fails. 

The main reason for this failure is that the probability of low 
mass clusters forming massive stars has values around 0.2 in the 
case of a 100 Mq zero age stellar cluster from our simulations. 
We note that this value must be understood statistically rather 
than interpreted as a value applicable to individual clusters. The 



correct interpretation of our findings is that 20% of 100 Mq zero 
age stellar clusters have stars massive enough to develop an H n 
region (i.e. they have at least one star more massive than 20 Mq) 
and 80% of zero age 100 Mq stellar clusters do not have massive 
stars and do not form any H n region at all. When the clusters are 
considered individually, this situation poses a severe problem, 
but when we observe a galaxy of high enough mass the mean 
value obtained by synthesis models for such a mass provides a 
good description of the system. 



However, a common misunderstanding of this result is that 
an individual cluster has 0.2 stars more massive of 20 Mq. Since 
the number of stars must be a natural number, it can be con- 
cluded erroneously that low mass clust ers do not form mas- 
sive s tars (the richness effect proposed bv lGarcfa Vargas & Dfazl 
1994), or, equivalently, that the upper mass limit of the IMF, 
m up , depends on the cluster mass, excluding the cases of 
clusters formed from mo lecular clouds less massive than m up 
dWeidner & Kr ouna 2006, and posterior works). This naive in- 
terpretation would imply that the ionizing flux (or, equivalently, 
the Ha luminosity) is not a good proxy of the star formation 
of galaxies, since an ad hoc cut-off in the upper mass limits 
of the IMF and the initial cluster mass fu nction (ICMF) is re- 
quired (c.f. iPflamm-Altenburg et~ai]l2009l) . In terms of g(H°), 
most low-mass clusters produce less ionizing flux than the mean 
obtained by synthesis models, but there are also a few low-mass 
clusters that produce a larger ionizing continuum than the mean 



11 



M. Villaverde, M. Cervino & V. Luridiana: Individual stars vs. stellar ensembles 



given by synthesis models, and the effect cancels out on average 
provided there are enough low mass clusters^. 

There are several observational results that are inconsistent 
with variations in the IMF upper mass limit, m up , as a func- 
tion of the cluster mass and confirm the probabilistic interpre- 
tation of the IMF (e . g. ICorbelli et al.l 120091 Uamet etalj|200l 
Mafz Apellani z et al.l 20071 as some examples). Even more, 
even if there were clusters with masses M c i us lower than m up 
(which implies a hard physical condition for the maximum stel- 
lar mass in the IMF realization), ;« up would be statistically 
preserved if the cluster system produced by the same molecu- 
lar cloud cont ains clusters with M c i us h igher than m up (size-of- 
sample effect. ISelman & MelnTckll20 08). 

With our photoionization simulations, we have illustrated the 
diversity of spectra that the scatter in r^ us implies and its im- 
plications in diagnostic diagrams. The large scatter shown in the 
diagnostic diagrams is caused only by the way in which we com- 
bine stars to obtain the ionizing continuum, which is compatible 
with a poor IMF sampling. They are not due to metallicity or age 
effects since the same abundances and ages have been used for 
all the simulations. When a comparison with theoretical grids of 
H ii regions based on the average ionizing continuum produced 
by synthesis models is performed, our theoretical clusters may 
be found to span a wide range in metallicity and/or age. 

4.1. The best ionizing continuum for low mass cluster 
modelling 

Our results have shown that obtaining the global properties of 
low mass clusters, such as their mass or age, is more difficult 
than expected, since the emission is dominated by individual 
stars. This situation turns out to be an advantage if we are in- 
terested in obtain a good description of the ionizing continuum 
of the cluster regardless of its evolutionary properties, since the 
properties of individual stars are better known than the statistical 
properties of clusters. 

In the particular case of the shape (parametrized as T'™*) and 
luminosity (parametrized as Q(H )) of the ionizing continum, 
the average obtained by synthesis models is not found be an op- 
timal representation of low mass clusters. At low masses, the 
cluster 7^ us and Q(H°) are strongly correlated and follow the 

T*£ ef - <2(H°) relation of individual stars. That is, the ionizing 
continua of low mass clusters (M c i us < 10 4 M0), including the 
case of unresolved ones, are statistically dominated by individ- 
ual stars. 

As a rule of thumb, we recommend the use of individual stars 
(with its corresponding stellar r* ff and Q(H°) values) rather 
than stellar ensembles (such as those implicit in synthesis mod- 
els) to model the ionizing flux of low mass clusters (M c \ us < 10 4 
Mq). However, this ionizing flux does not provide information 
about the global properties of the cluster (especially the clus- 
ter mass or age). This rule can also be restated in the following 
terms: ionizing low mass clusters (including unresolved ones) 
are more suitable for massive stellar atmosphere studies (ioniz- 
ing spectra studies) than for evolutionary ones. 

We note that the relation that links r| us to Q(He°)/Q(H°) 
has been derived by assuming ZAMS models and a specific 
set of atmosphere models. A different choice would lead to 



4 Although it is the case for properties such as Q(H°), Ha, and stel- 
lar luminosities whose average values scale linearly with the cluster 
mass, it would not be the case for properties that scal e in a non-linear 
way, such as the intensity of collisional emission lines (Villaverde et al. 
2009). 
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Fig. 14. T* vs. g(He°)/g(H°) for individual stars assumed 



in this work based on no n-LTE blanketed CoStar models 
(Scha erer&deKoterl 1 19971) (soli d line X iMas-Hesse & Kunfhl 
d 19911) relations based on iMihalasI (1 19721) non-LTE unblanketed 
models (red dotted line) for stars with luminosity class V a nd I 
(hexagons and triangles respe ctively), Martins et al. (2005|) re- 
lations based on CMFGEN dHillier & Milled 1 19981) non-LTE 
spherical extended line-blan keted mod e ls (bl ue short-dashed 
line) for V a nd I stars, and Smith et al] d2002l) relations based 
in WMbasic (Pauldrach et al .1 120011) computations (orange long- 
dashed line). 



a different reference scale: generally spe aking, any reference 
scale is atmosphere - model dependent (see | Morisset et al.| [2004, 
iMartins et alJ l2005l ISimon-Dfaz & Sta sihska 2008, as exam- 
ples). However, any reference scale would be qualitatively simi- 
lar and lead to similar results that depend on only two results: (a) 
the IMF in low mass clusters is poorly sampled (so the number of 
ionizing stars is intrinsically low, and it is just one in low mass 
ionizing clusters) and (b) the relation between the stellar mass 
and its ionizing properties is non-degenerate (in such a way that 
sampling effects in the mass translate directly into the ionizing 
properties). As long as (b) holds for a particular calibration, the 
outcome will be the same. This is generally true, as illustrated 
by Fig. M 

Although our reference scale has been compiled with par- 
ticular stellar models, we recall that T^ s is not the T* B of any 
particular star, but rather a convention al measure of the cluster 
Zanstra-like temperature dZanstral 19271) . That is, the scale can be 
used as a conventional one and applied to evolved clusters, even 
if in these the most luminous stars are not ZAMS stars and the 
stars they contain do not necessarily fit into the reference scale 
(such as Wolf-Rayet ones). But even if we wish to account for 
evolution, Fig. [14] shows that explicit consideration of different 
luminosity classes would only add some dispersion to the scale. 
In sum, our scale is conventional but not arbitrary. 

In practice, there may be additional difficulties to consider, 
such as the presence of more than one ionizing star in the unre- 
solved cluster. However, even in these cases the correct approach 
is not using the mean value from synthesis models, but rather add 
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the contribution of several ionizing stars, in such a way that the 
t:„ and Q(H°) of the mix reproduces the observed values (e.g. 
Jamet et al. 200j|). 



5. Conclusions 

We have studied the distribution of stellar cluster properties for 
a given age and metallicity and the influence that the sampling 
effects of the IMF may have on it. We have also explored the 
implications for Hn region spectra and the estimation of clus- 
ter and H n region properties obtained with grids and diagnostic 
diagrams. 

For this task, we have performed 2x 10 6 Monte Carlo simula- 
tions of zero age clusters in the cluster mass range 20 - 10 6 Mq 
to estimate the distribution of the cluster properties. Focusing 
on clusters with masses around 100 Mq (the most numerous 
ones), we conclude that only 20% of these clusters can gener- 
ate an H ii region. We also show that the shape of the average 
spectra produced by synthesis models only represents 8% of the 
possibles cases (or 33% of clusters that can develop an H n re- 
gion). However, those clusters have a higher ionizing flux lumi- 
nosity than the one produced by the average value obtained by 
synthesis models, so the average values obtained from synthesis 
models is correct statistically when all possible cases are consid- 
ered, but not for individual clusters. 

We have shown that variations in the upper mass limit of the 
IMF with the cluster mass mimics strong sampling effects for 
some (but not all) stellar clusters. However, a correct represen- 
tation of low mass star clusters is not compatible with an ad hoc 
dependence of the upper mass limit of the IMF on the cluster 
mass since it does not cover all possible cases. 

We have also found that there is a strong correlation between 
the cluster ionizing flux (parametrized as Q(H )) and the shape 
of the ionizing continuum produced by the cluster. This allows us 
to describe the ionizing flux of a low mass cluster (M c i us < 10 4 
Mq) as that produced by a individual star, instead of using an 
ensemble of stars. Thus H n regions with ionizing stellar cluster 
masses lower than or about 10 4 Mq are suitable for direct studies 
of individual hot-star atmosphere studies even when individual 
stars are not fully resolved in the cluster. 

Using photoionization models, we have shown that it is not 
useful to use grids of photoionization models to obtain global 
properties of low mass clusters, since they produce erroneous 
results. 

Finally, we have illustrated how IMF sampling effects work 
and the reasons why the modelling by using individual stars of 
H ii regions produced by low mass clusters is more correct than 
the use of the average parameters obtained by synthesis models 
such as the cluster 7^ us : sampling effects are a trade-off between 
global properties of the cluster and particular properties of the 
stars. When the observables are determined by only a few stars, 
global properties are difficult to obtain, and when the observ- 
ables are defined by the ensemble as a whole, the information 
about particular stars are diluted by the ensemble. 
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